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ABSTRACT 



In a previous paper we have written down equations describing steady-state, 
optically thin, advection-dominated accretion onto a Kerr black hole ( |Gammie 
fc Popham (199T)| , hereafter Paper I). In this paper we survey the numerical 
solutions to these equations. We find that the temperature and density of the 
gas in the inner part of the accretion flow depend strongly on the black hole 
spin parameter a. The rate of angular momentum accretion is also shown to 
depend on a; for a greater than an equilibrium spin parameter a^q the black 
hole is de-spun by the accretion flow. We also investigate the dependence of 
the flow on the angular momentum transport efficiency a, the advected fraction 
of the dissipated energy /, and the adiabatic index 7. We find solutions for 
-1 < a < 1, 10"^ < « < 0.44, 0.01 < / < 1, and 4/3 < 7 < 5/3. For low values 
of a and / the inner part of the fiow exhibits a pressure maximum and appears 
similar to equilibrium thick disk solutions. 



1. Introduction 

In advection-dominated accretion fiows (ADAFs) the accreting gas fiows inward much 
more rapidly than it can cool. The energy released by accretion goes into heating the gas. 
This is in marked contrast to the usual thin accretion disk, where the radial velocity is 
small and the accretion energy is efficiently radiated away. If a black hole accretes via an 
ADAF much of the accretion energy can be carried across the horizon with the heated 
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gas, reducing the luminosity well below that of a comparable thin disk. Paper I briefly 
summarizes the development of advection-dominated disk theory. Over the past few years, 
this theory has had notable success in reproducing the observed spectra of black hole 
candidate systems. For a review of the theory and applications of advection-dominated disk 
theory, see [Narayan (1997) . 



Early models of advection-dominated flows around black holes (|Chakrabarti (1996) 



Nakamura et a/.(1996)| , [Nakamura et ffl/.(1997)| , |Chen, Abramowicz, fc Lasota (1997) 



Narayan, Kato, fc Honma (1997)| ) did not include a proper treatment of relativistic effects. 



but instead used an approximate pseudo-Newtonian potential due to [Paczyhski fc Wiita 



(1980)| . Close to the black hole event horizon, the gas temperatures and velocities can 
become extremely high. This hot, rapidly rotating and infalling gas should produce 
important observable effects in the high-energy spectra of black hole candidates. Relativistic 
effects not accounted for by the pseudo-Newtonian potential are dominant in this innermost 
region. The character of the inner portion of the disk also depends strongly on the black 
hole spin a, which is not included in the pseudo- Newtonian treatment. For these reasons, it 
seems clear that a fully relativistic model is needed. 

Recently, advection-dominated disk models have begun to include the effects of general 
relativity, allowing a more accurate examination of the flow close to the event horizon. 
Both [Abramowicz et al. (1997] (AGCL) and [Peitz fc Appl (1997)| (PA) have presented 



disk solutions in the Kerr metric. In Paper I, we wrote down a set of disk equations in the 
Kerr metric, and showed a few example solutions. Our relativistic disk equations differ in a 
number of respects from those of ACGL and PA. These differences are described in detail 
in Paper I. The main differences with ACGL are that we use a causal stress and include 
the relativistic enthalpy. PA also include these effects, although they use a more simplifled 
prescription to enforce causality. Unlike PA, who use a polytropic equation of state, we 
solve the energy equation assuming that a constant fraction of the dissipated energy gets 
advected with the gas. Also, unlike ACGL and PA, we use the height prescription of 
Abramowicz, Lanza, fc Percival (1997) . 



In this paper, we examine the structure of our numerical solutions for steady state 
flow onto a Kerr black hole in detail. In particular, we explore the effects of changing the 
dimensionless black hole spin a, the viscosity parameter a, the advected fraction of the 
dissipated energy /, and the adiabatic index 70 (ACGL have shown solutions for three 
values of a, while PA presented a more extensive set of solutions for various values of a and 
a). We show that there is an equilibrium spin rate for black holes accreting via ADAFs, 
above which the black hole will be spun down by accretion. We also show that there is a 
pressure maximum in the vicinity of the last stable orbit for low a solutions (cf. NKH). 
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Finally we show that the Bernoulli parameter can be positive in some regions of the flow, 
which at least suggests the possibility of a pressure-driven outflow. 

To familiarize readers with our notation, we summarize our relativistic disk equations 
and method of solution in §2. In §3 we examine the dependence of our disk solutions on 
a, a, /, and 70. §4 demonstrates that there is an equilibrium spin rate for black holes 
accreting from ADAFs. We discuss the implications of our results in §5. 



2. Disk Equations and Solution Method 

2.1. Summary of Disk Equations 

We use a relativistic version of the slim disk equations first introduced by |Paczyhski| 
fc Bisnovatyi-Kogan (198T]] and [Muchotrzeb fc Paczyhski (1982)| to study the structure 



of advection-dominated accretion disks. The equations are presented in detail in Paper 
I, but for the reader's convenience we summarize them here. The units are such that 
G = M = c = 1, where M is the mass of the black hole. 

We begin with the continuity equation: 

Arcr'pHeV (j— ^) = -M. (1) 

Here r is the Boyer-Lindquist radius, p is the rest mass density, = H/r is the relative 
disk thickness, V is the radial velocity measured in a corotating frame, T> = 1 — 2/r + jr^ 
is a relativistic correction factor, and M is the rest mass accretion rate. 

The gas energy equation is 

v.f^V'7|!if (2) 



) \ dT dr p dr ^ 

Here u is the internal energy per unit proper volume, T is a dimensionless temperature such 
that the pressure P = pT, $ is the dissipation function, and / < 1 is a parameter used to 
mock up the effects of cooling, if any. The dissipation function $ is given in Paper I (eq. 
(68) and Appendix A). We also take 

which is a good approximation to the exact relativistic equation of state for an ideal 
relativistic Boltzmann gas (see [Chandrasekhar (1939)| ) when 70 = 5/3. Because the 
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accretion flow is assumed to consist of a mixture of magnetic fields and ionized plasma, 70 
can be smaller than 5/3. 



The radial momentum equation is 



(4) 



1 — dr pr] dr 

where 

The fr term combines the effects of gravity and rotation, where ^ = 1 + a^/r^ + 2a^/r^ 
and 7^ = 1 + — V'^)/{r'^A), ^ = u'^/u* is the angular velocity, and fi± = ±(r^/^ ± a)^^. 
The radial acceleration, on the left-hand side, is given by the difference between fr and the 
pressure gradient force, where i] is the relativistic enthalpy r] = {p + p + u) / p. 

The angular momentum conservation equation is 

MIt] - ATiHerH''^ = Mj. (6) 

Here / is the specific angular momentum of the accreting gas, j = const, is the angular 
momentum accretion rate per unit rest mass accreted. The remaining term gives the viscous 
angular momentum transport rate, where is the viscous stress. Notice that we solve for j 
as an eigenvalue, as do NKH and ACGL. PA prespecify j (which they label Lq) and iterate 
in j when they want to find a solution with a specified value of I at the outer edge of the 
disk. |Chakrabarti (1996]| also prespecifies j. 



The calculation of is rather lengthy, and the reader is referred to §4 of Paper 1 
for a full mathematical discussion. To briefly summarize the relevant physics, we begin 
by calculating the shear rate a in the local rest frame of the accreting gas. This gives a 
complicated expression which is closely approximated by the thin disk shear rate derived 
by [Novikov fc 'I'horne (1973)| , crthin = {I / 2)A'y^rdQ / dr (see the Appendix of Paper I for a 
discussion). Next we specify the relation between the turbulent shear stress and the shear 
rate (rate of strain). The simplest prescription is the usual Navier-Stokes form in which 
stress is simply proportional to rate of strain. As is well known, the resulting system of 
equations is acausal. To preserve causality, we use a relativistic version of a prescription 
that originated with [Maxwell (1867)| and Cattaneo (1948, 1958) and was discussed in the 



context of accretion disk theory by [Papaloizou fc Szuszkiewicz (1994]| ([Narayan, Loeb, 



Kumar (1994)| give a mathematically equivalent formulation). This modification limits the 
propagation of viscous effects to a finite speed Ci,, which we take to be Ci, = -JaCg (our 
solutions are not sensitive to the precise value of Cy). The viscous stress is then obtained 
by transforming back from the local rest frame to the Boyer-Lindquist frame. 
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The equation of vertical mechanical equilibrium is 



Hi = (7) 



where z/^ is an effective vertical frequency. We adopt the expression derived by Abramowicz, 
Lanza, fc Percival (1997)| for z/^: 



where S = —Ut is the "energy at infinity" , which is conserved along geodesies. 
Finally, the inward flux of mass-energy is 

E = A'KHeA-{p + u + p)£vJ- + tt''l (9) 

where t is the viscous stress tensor. We do not use this equation to flnd our solution. But 
when / = 1, is constant with r, and E ~ M, because T <^ 1 at large r. In this case 
equation is a check on our numerical accuracy. It is generally satisfled to better than 
one part in 10"'^. For / < 1, however, E is not conserved, since some of the rest-mass is 
radiated away. Recall that E at the horizon, and not M, is the true rate of change of the 
mass of the black hole. 



2.2. Critical Points and Boundary Conditions 

We set the outer edge of the ADAF at 10^ Schwarzschild radii, or r = 2 x 10"^ GM/c^, 
and the inner edge just outside the event horizon at r = (1 + yl^^a?) GMj(?^ where P = 0. 

At the outer edge we impose two boundary conditions: and must equal their 
values in the self-similar advection-dominated solution of [Narayan fc Yi (1994)| . We have 



also obtained solutions with the thin disk values for f2 and Cs at the outer edge, but 
our solutions adjust to the self-similar proflle within a short distance of the boundary. 
NKH found a similar result, since their equations are essentially equivalent to ours in the 
non-relativistic limit. Thus the self-similar solution is the "natural" state of the flow far 
from the black hole, and is an appropriate outer boundary condition. 

Two other conditions on the flow are provided by the requirement that the flow pass 
smoothly through two critical points. The flrst is the sonic point r^, where \y\ ~ c^. The 
second is the "viscous point" r„ associated with the flnite propagation speed of viscous 
effects, where \y\ ~ Cy. Associated with each critical point are two conditions that must 
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be satisfied for a smooth fiow, as well as one degree of freedom, the location of the critical 
point itself. 

The final boundary condition normalizes the density (the density appears in the basic 
equations only in the form dlnp/dr). For simplicity we set the normalization constant so 
that M — 1. We now have all the boundary conditions required to solve the four first-order 
ordinary differential equations for V,l, p, and T, and to find the eigenvalue j. In particular, 
consistent with causality, no boundary conditions are applied at the event horizon. 

2.3. Method of Solution 

We solve the system of equations listed above using a relaxation method. We solve for 
the structure of the flow in three radial regions: the subsonic, "sub-viscous" outer zone, 
where r > r„, the subsonic, "super- viscous" middle zone, where < r < r^,, and the 
supersonic, "super- viscous" inner zone, where r < Vg. We first solve the outer zone for 
a specified value of r^, and then solve the middle zone for the same r^. We compare the 
outer and middle zone solutions to sec whether they match up at r^; if not, we change r„ 
and repeat the procedure until the outer and middle solutions match up. Then, using the 
variable values from the middle solution at as boundary conditions, we solve for the inner 
zone. 

As is well known, our independent variable, the Boyer-Lindquist r, is ill-behaved near 
the horizon when \a\ —>■ 1. This might raise some concerns about the numerical accuracy 
of our solutions. We have taken several steps to control for this, however. First, we use 
a fixed number of grid points in each radial region. As a — > 1, the sonic radius moves 
sharply inward, and this greatly increases our numerical resolution close to the horizon. 
Second, we have convergence-tested our solutions by increasing the number of grid points. 
No significant changes were found at increased resolution. Finally, notice that conditions 
close to the horizon do not influence the solutions at larger radius because our fundamental 
equations arc causal: the solution in the inner region is obtained by integrating inward from 
the sonic point. In fact, many of the important parameters of the flow, such as the angular 
momentum accretion rate j, are fixed by conditions at and outside the sonic point. 

3. A Survey of Solutions 

Here we describe the changes in our solutions that result from changes in four basic 
parameters: the black hole spin a, the viscosity parameter a, the advected fraction of the 
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dissipated energy /, and the adiabatic index 7. We begin from a base solution with a — 0, 

a — 0.1, f — 1, and 7 = 1.4444 (which corresponds to equipartition between gas and 
magnetic pressure). This sohition was described in detail in Paper I, and it appears in each 
of the sets of solutions described below. Recall that physical units may be recovered as 
follows: radial velocity is Vc, angular momentum is IGM/c, density is pMG/c" (since the 
mass accretion rate is set to 1), and temperature is Trhc^/k, where rh is the mean molecular 
weight and k is Boltzmann's constant. 

It is worth pointing out immediately that the marginally stable orbit plays no special 
role in many of our solutions. While the sonic point often lies close to the radius of the 
marginally stable orbit, there is no abrupt change in the character of the flow there. This is 
because pressure gradients play an important role in the radial structure of ADAFs; unlike 
thin disks, ADAFs are not in centrifugal balance. 

3.1. Black Hole Spin a 

Wc begin by examining the dependence of the accretion flow on the black hole spin 
parameter a. We solve for values in the range — 1 < a < 1, allowing for the possibility of a 
counterrotating hole. The horizon lies at r = 1 + (1 — a^)^/^, i.e. at r = 2 for a = and at 
r = 1 for a = ±1. The resulting profiles of p, V, fi, Z, T, and H/R are shown in Figure 1 
for solutions with a = -0.999, -0.9, -0.5, 0, 0.5, 0.9, and 0.999. 

Notice that changes in a have little effect far from the hole. Also, notice that in 
general the positive-a solutions differ from the a — solution more than the corresponding 
negative-a solutions. This is especially dramatic in p and T: the a — 0.999 solution reaches 
very high density and temperature near the horizon, while the a — —0.999 solution has p 
and T at the horizon similar to the a — solution. 

The positive and negative-a solutions also show very different V profiles close to 
the hole. The positive-a solutions accelerate rapidly and plunge into the hole, while the 
negative-a solutions accelerate rapidly farther from the hole and only gradually near the 
hole, coasting across the horizon. 

Figure 1 also shows the positions of the sonic and viscous points for each solution. The 
sonic and viscous points move to smaller r as a increases, as expected since the "plunge" 
toward the hole occurs at smaller radii. At a = 0, the sonic point is at r ~ 6.41, close to 
the last stable orbit. At a = 0.999, it has moved in to r ~ 1.94, while for a = —0.999 it 
moves out to r ~ 9.20. The viscous point always lies much farther from the hole; in general 
r„ ~ 4 — Bts- This is due to our setting = ac^, so that for a = 0.1 we have Cj, ~ O.Sc^. 
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The value of at the horizon should be equal to a; = 2a/ {r^ + (j?r + 2a^) (denoted by 
a dashed line in Fig. 1) due to the effects of frame- dragging. This results in the symmetric 
distribution of inner values between -0.5 and 0.5 seen in Fig. 1. For all but the high-a 
solutions, f2 shows a boundary layer-like profile, reaching a peak and then dropping rapidly 
to uj at the horizon. Had we plotted O = O — a;, the profiles for all values of a would 
resemble the ^2 profile for a = 0. 

The specific angular momentum Iri drops steadily as a increases. Notice that Iri includes 
the angular momentum associated with both the mass density and the energy density of 
the gas. The inner value of Iri decreases from over 3 to less than 2 as a goes from -0.999 
to 0.999. The actual rate of angular momentum accretion by the hole per unit rest mass 
accreted is j. Values for j are generally close to the horizon values of Ir]. This means that 
the rate of viscous angular momentum transfer is very small near the horizon; however, it 
is not zero. 

The relative disk thickness Hjr also varies dramatically with a. Notice that for all 
values of a, Hjr is smaller near the horizon, then increases to ~ 1 at large radii. Negative-a 
flows are much thinner than positive-a flows. 

3.2. Viscosity Parameter a 

Figure 2 shows a set of solutions with a = 0.001, 0.003, 0.01, 0.03, 0.1, 0.3, aU with 
a = 0, / = 1, and 7 = 1.4444. Both the sonic and viscous points move steadily outward as 
a increases, since a larger a gives a larger viscosity coefficient for given Cg and H, producing 
larger radial velocities V. This removes more angular momentum and dissipates more 
energy, resulting in smaller values of Irj and higher temperatures for larger a. The higher 
temperature gives a larger relative disk thickness H/r, and this combined with the larger 
radial velocity makes the density decrease as a increases. 

The density and temperature proflles show an interesting effect at low a: the 
development of a local maximum in density and temperature just outside the last stable 
orbit. This effect was also noted by NKH and CAL. The maximum becomes more 
pronoimced as a decreases. It is faintly seen as an inflection in the density profile at 
a = 0.01, but by a = 0.001 it produces a very large, wide peak with a density 50% larger 
than the density at the horizon. The maximum in the temperature profile is less dramatic, 
but the combined effect of the two is to produce a pressure maximum which affects the 
dynamics of the flow. On the inner side of this pressure maximum, the pressure gradient 
force points inward. This inward force is balanced in part by more rapid rotation: Q 
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increases as a decreases. In fact the a — 0.001 solution even has a super- "Keplerian" 
rotation in a small region extending from r ~ 5 — 8. 

3.3. Advected Fraction / 

The parameter / was introduced to mock up the effects of cooling. The case of greatest 
interest for modehng ADAFs is when cooling is unimportant, so / ~ 1. We can reduce / 
so that it is much less than 1, however, and in this case we should approximately recover a 
thin disk, albeit with a radial temperature structure that is not consistent with any realistic 
cooling function. 

Figure 3 shows solutions for f — 1, 0.3, 0.1, 0.03, 0.01. All of the solutions have a = 0, 
a = 0.1, and 7 = 1.4444. The dramatically lower temperatures in low-/ solutions reduce Cg 
and H and thus reduce the viscous stress. As a result, this sequence of solutions strongly 
resembles the a sequence shown in Fig. 2. For small values of /, as for small values of a, a 
maximum develops in density, temperature, and pressure in the inner disk for small values 
of /, resulting in super-Keplerian values of O. The relative disk height drops as / decreases, 
and the removal of angular momentum is less efficient. The major difference between the 
two sequences is that low values of / produce much lower temperatures and higher densities 
than low values of a. 

The low-/ solutions are indeed close to a thin disk. In particular, the sonic point 
occurs close to the last stable circular orbit (this is true for solutions with a variety of a at 
low /). Also, the energy-at- infinity £ and angular momentum I are close to their circular 
orbit values. Of course, this is expected due to the low temperatures of these models, so 
that radial pressure gradients play a negligible role. Nevertheless, it is good to see this 
expectation confirmed. 

3.4. Adiabatic Index 7 

The adiabatic index is not fixed by the temperature because some (unknown) fraction 
of the total pressure is contributed by the magnetic field, while the rest is contributed by 
gas pressure. Figure 4 shows a sequence of solutions for a range of values of 7 = 1.3333, 
1.4444, 1.55, 1.66. These correspond to the fraction of the pressure contributed by the gas 
pressure /? = 0, 0.5, 0.8, 0.99, respectively. All solutions have a — Q,a — 0.1, and f — I. 

The high-7 solutions get much hotter than the low-7 ones, increasing the effective 
viscosity coefficient. This produces higher velocities, larger relative disk heights, and 
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removes more angular momentum. As 7 approaches 5/3, these effects become extreme. For 
example, the 7 = 1.66 solution has unphysically large values of H/R. 



4. Equilibrium Black Hole Spin Rates for Advection-Dominated Accretion 

The black hole accretes both mass and angular momentum, so the dimensionless black 
hole spin parameter a = Jc/GM"^ (J is the black hole angular momentum) changes with 
time. In dimensionless form, 

^ = Mj - 2aE, (10) 

where E is the inward flux of mass-energy. Thus if j = 2aE/M, then da/dt vanishes and 
the black hole has reached an equilibrium spin Ogg. The value of Oeg, if it exists, will depend 
on the input parameters a, a, /, and 7. Notice that when / = 1, as in most of the solutions 
presented here, E = const. ~ M, and Ogg is reached when j = 2a. 

The equilibrium spin is of considerable astrophysical interest since the observable 
properties of accretion flows change sharply in the neighborhood of a = 1. In an influential 



paper, Bardeen (1970) showed that a black hole accreting from a thin disk reaches a = 1 



in a finite time. [Thorne (1974)| pointed out that if one takes account of the preferential 



accretion of angular momentum by the hole, then a^q = 0.998. [Abramowicz fc Lasota| 



1980)1 suggested that even larger values of Ueq might still be obtained from a thick disk, 
based on numerical calculations that showed that the inner edge of these disks was located 
somewhere between the last stable orbit and the last bound orbit. 

In order to calculate Ogg, we have obtained sequences of solutions with increasing a. As 
a increases, j decreases, and when j = 2aE/M we have found a^q. In Figure 5 above we 
show the dependence of j on a for the six values of a = 0.001, 0.003, 0.01, 0.03, 0.1, 0.3, for 
which a = solutions were shown in Fig. 1. In all cases, j decreases as a increases, and 
decreases faster as a approaches 1, particularly in the low-a solutions. Here / = 1, so 
is reached when j = 2a, as denoted by the dotted line. Notice that aeq is quite far from 
a = 1 for high-a flows: for a = 0.3, we have a^q — 0.8. On the other hand, for low-a flows, 
is very close to 1: for a = 0.001, we have aeq = 0.999965. The variation of aeq with a 
is shown in Figure 6. For the sequence of a = 0.001, 0.003, 0.01, 0.03, 0.1, 0.3, we find 
aeq = 0.999965, 0.99972, 0.9975, 0.985, 0.930, 0.806. For a < 0.03 we find the approximate 
relation 1 — aeq — lOa^'^. 

We also obtained sequences in a for values of / between 0.01 and 1, using 
a = 0.1,7 = 1.4444. For / < 1 we no longer have E ^ M, so the full form of 
equation (0) is required. We find aeq = 0.930, 0.988, 0.9984, 0.99983, 0.999979, for 
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/ = 1,0.3,0.1,0.03,0.01, respectively, as shown in Figure 6. Here again we find an 
approximate power-law relationship 1 — Qeq — 0.1/^'^. 

These results have some interesting implications. If most of the mass of a black hole is 
accreted through an ADAF, then the hole cannot approach a = 1. All advection-dominated 
models of observed systems that have appeared in the literature have high values of a ~ 0.3, 
since solutions with low a have extremely low luminosities. According to the results given 
above, the black holes in these systems will not be able to spin up past a ~ 0.8. The most 
dramatic effects of the black hole spin on the flow, which appear when a nears unity, will 
not be seen in such systems. 

Of course, the accretion of photon angular momentum is not included in our calculation, 
but this is justified where 1 — / is small. In our flows, small values of / correspond to thin 
disks, which radiate a fraction 1 — / of the energy dissipated in the disk. The limiting spin 
aeq also depends on /, and we find that flows with / > 0.1 will have a^q < 0.998, so they 
should stop spinning up before reaching the photon-transport limit. Flows with / < 0.1 
have tteq > 0.998, so here the photon-transport limit is relevant. 

5. Discussion 

5.1. Comparison to Earlier Solutions 

Our disk equations differ in a number of respects from those used by other authors, 
as detailed in Paper I and above. Also, the solutions described above cover a wider range 
of parameter space than previous studies. We have found solutions over the ranges of 
a = -0.99999 to 0.99999, a = 0.0001 - 0.44, / = 0.01 - 1, and 7 = 1.3333 - 1.66. 

Our sequence of solutions for a = 0.001 — 0.3 resembles the solutions of NKH in many 
respects. We find sonic radii ranging from ~ 4.28 — 9.41 as a increases, similar to the 
range = 4.20 — 10.63 found by NKH. Our solutions also develop pressure maxima at 
small values of a, as discussed below. Since NKH used the pseudo- Newtonian potential, 
their solutions have some unphysical properties near the horizon, which are avoided by our 
use of the full relativistic equations. 

Our relativistic solutions with a = produce emission spectra which are qualitatively 
similar to those calculated from the NKH solutions. The densities tend to be somewhat 
higher near the horizon in our solutions due to the inclusion of relativistic effects, and this 
increases the luminosity of the flow substantially. This was shown by |Narayan et ffl/.(1998)| , 
who calculated the emission spectrum of Sgr A* using the relativistic dynamical solutions 
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described in this paper. They were thus able to make a direct comparison between the 
spectrum computed from our solution and that computed from a pseudo-Newtonian solution 
as described by NKH. This showed that the inclusion of relativistic dynamics increased 
the emission by a factor of 10-100 in the infrared. Note that relativistic photon transport 
effects are not included in these spectra, apart from including the gravitational redshift. 
Also, the effects of using rotating black hole solutions have yet to be explored; however, it 
is clear that the substantial increases in temperature and density near the horizon in these 
solutions will have a dramatic effect on the emission spectrum. 

ACGL presented solutions with a = 0, 0.5, and 0.99 for a = 0.1, which included 
bremsstrahlung cooling. The shape of the /(r) curves is similar to ours; however, the 
solutions differ in some other respects. First, the ACGL solutions have larger values of I 
at the inner edge of the flow. They find Zj„ ~ 3.2, 2.6, 1.7 for a — 0, 0.5, 0.99, respectively, 
whereas our solutions have Zj„ = 2.14, 1.79,0.93 for these same values of a. One reason 
for this is that their l{r) curves only extend in to the Schwarzschild radius at r = 2, even 
for the a = 0.5 and a = 0.99 solutions, where the horizon sits at smaller radii. Also, in 
ACGL's scheme, the angular momentum eigenvalue which we call j is equal to 1^ (or in 
their notation, C = Cq), whereas in our solutions, j can differ substantially from /j„ due 
to viscous torques. Our values of j for these solutions are 2.62, 2.31, and 1.73, which are 
somewhat closer to the ACGL values. Another difference between our solutions and those 
of ACGL is that their solutions show a maximum in log Cg and inflections in log P close to 
the inner edge, whereas our solutions have log T and log P increasing smoothly all the way 
in to the horizon. 

PA have presented the most extensive survey of solutions to date. They find two types 
of viscous disk solutions: "type I" solutions at low values oi a {a — 0.001 — 0.045 in the 
solutions shown) and "type II" solutions at high a {a > 0.3 in the solutions shown). One 
major difference between PA's scheme and ours is that they specify the angular momentum 
eigenvalue j (which they call Lq) while we solve for it. Their type I solutions with a ranging 
from 0.01 to 0.045 all have the same value of Lq, and show great variation in their angular 
momentum profiles at large r, whereas our solutions tend to have rather similar values 
of / at large r. PA's type 1 solutions also have a maximum in the sound speed outside 
the sonic point, and the sound speed then decreases down to the horizon. This may be 
a consequence of their polytropic equation of state, which requires that the sound speed 
must decrease inward if the radial velocity increases inward more rapidly than ~ r~^. Their 
lowest-a solutions a,t a — 0.001 — 0.008 have a super-Keplerian region in the inner part of 
the disk. In the type II solutions shown by PA, the sonic point occurs far out in the flow at 
r ~ 30 — 60. This is a much larger sonic radius than in our high-a solutions: our a = 0.3 
solution has = 9.41. According to PA, their type II solutions have the sonic point located 
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at an outer critical point rather than an inner critical point as the Type I solutions do. 
Unlike PA, who found two types of solutions, and [(Jhakrabarti (1996)|, who has found 



solutions where the flow goes through radial shocks, we find smooth solutions for the entire 
range of parameter values, and we see no evidence for sudden transitions between different 
types of solutions. In order to look for alternate solutions, we varied r„ and solved for the 
outer and middle sections of our solutions to see whether these sections would match up 
for an alternate value of r„. Despite looking for alternate solutions with several choices of 
parameters, including low values of a, we found no additional solutions. 

It is worth commenting in some detail on why our solutions are shock-free while those 



of [Chakrabarti (1996)| are not. The origin of this difference lies in how the boundary 
conditions are treated. The difference is most easily explained by analogy with the problem 
of spherical (Bondi-Hoyle) accretion. To solve the Bondi problem one specifies the density 
and temperature, but not the radial velocity (nor, equivalently, the accretion rate), at large 
distance from the accreting object. The velocity at large radius is adjusted until the flow 
passes smoothly through the sonic point. This velocity, or the accretion rate, is thus an 
eigenvalue of the problem and must be solved for self-consistently. This approach has been 
validated by numerically solving a realistic initial value problem and showing that it settles 
down to the Bondi solution. If one were to specify the radial velocity or accretion rate the 
flow would not generally pass smoothly through the sonic point; rather, it would shock. 

Our treatment is analogous to the standard treatment of the Bondi problem, except 
that in our case the mass accretion rate is specified, while the angular momentum accretion 
rate j is the eigenvalue. It is adjusted so that the flow passes smoothly through the sonic 
point. Chakrabarti's approach, on the other hand, is analogous to prespecifying the mass 
accretion rate in the Bondi problem: he fixes j. As a result the flow does not generally 
pass smoothly through the sonic point and his solutions contain one or more shocks. Our 
intuition, and the analogy with the Bondi problem, suggest that this is not the correct 
approach. The issue can only be settled conclusively, however, by solving a realistic initial 
value problem and showing that it converges to one solution or the other. This has not yet 
been done. 



5.2. Local Pressure Maxima 

Our low-a solutions have a maximum in pressure, density, and temperature in the 
inner disk. Figure 7 shows density contours in the inner parts of the flow for various values 
of a and a. The contour plots are in the x = rsm{9),z = r cos{9) plane, and assume that 
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p oc exp(— 0^/(2iJg)). These solutions resemble the low-a solutions of NKH. They are 
also similar to the thick disk models developed by a number of workers (see [l''ishbone fc 



Moncriet [WlW% [Faczyhski fc Wiita (iMJJl , |Kees et ai. (1982)| ). In particular they have 



a pressure maximum, approximately constant specific angular momentum, and a region 
of super-Keplerian rotation where pressure support is important. Generally the pressure 
maximum lies close to the last stable orbit, and the sonic point lies close to the marginally 
bound orbit. Our approximations are most accurate for these low-a solutions in that H/R 
is relatively small (see Figure 2) and the low viscosity implies a disk with relatively low 
turbulent velocities. 

It is not entirely clear that the low-a models are relevant, however. A lower limit to the 
efficiency of angular momentum transport is set by the existence of the global hydrodynamic 
instability of [Papaloizou fc Pringle (1984)|. Simulations of the nonlinear development of the 



magnetorotational instability (Palbus, Gammie, fc Hawley (1995) ) suggest an even larger 



lower limit of order a ~ 0.01. For a as large as this a pressure maximum does not develop, 
and the pressure increases smoothly down to the event horizon. 

A new result of this study is the similarity of low-/ solutions to low-a solutions. Our 
/ = 0.01 and / = 0.03 solutions show pronounced density and pressure maxima in the 
same region of the flow, just outside the last stable orbit. They also have super-Keplerian 
rotation in this region. However, the low-/ solutions have the sonic point at ~ 5.72 for 
/ = 0.03 and at ~ 5.85 for / = 0.01, closer to the last stable orbit than the marginally 
bound orbit. 



5.3. Bernoulli Parameter 

The Bernoulli parameter Be measures the sum of the kinetic energy, potential 
energy, and enthalpy of the gas. Narayan & Yi (1994, 1995) pointed out that self-similar 
advection-dominated flows have a positive Bernoulli parameter for / > 1/3. The positivity 
of Be suggests the possibility of a pressure-driven outflow; Be is conserved for adiabatic, 
inviscid flows, so that gas with Be > could flow outward adiabatically and still have 
positive kinetic energy at large radius. 

In its relativistic form. Be = rjS — 1, where S = —Uf. Figure 8 shows Be for the 
solutions shown in Figs. 1-4, illustrating the variation of Be with radius and with a, a, /, 
and 7. Solutions with small values of a reach only small positive values of Be. In these 
solutions. Be tends to peak at r ~ 6 — 20, and reaches peak values of Be < 0.01 (see Fig. 
8b, d). Solutions with a approaching unity can have substantially larger values of Be which 
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peak at the horizon (Fig. 8a). Solutions with / < 1, which radiate away a fraction 1 — / of 
the dissipated energy, have Be < 0. At large radii, the solution with / = 0.3 has Be only 
slightly negative, as expected from the self-similar result that ADAFs with / > 1/3 have 
Be > 0. 

Variations in the Bernoulli parameter are directly related to the radial viscous energy 
flux, since the total radial flux of mass energy is conserved. We have 

ip j-r 

Be^^-1 + -^ (11) 

where t^j, is the viscous stress tensor. Recall that E at the event horizon is the actual rate 
of change of the black hole mass; furthermore, E = const. ~ M when / = 1. One can show 
that, if S is the shear stress measured in the local rest frame of the fluid, then 

tl = -rVu^S. (12) 



Since u"^ ~ at the horizon, t[ is finite at the horizon. Figure 9 shows the run of —t^, 
the outward viscous energy flux, for solutions with several values of a but otherwise with 
the standard parameters. 

The existence of a finite outward angular momentum and energy fiux at the horizon 
suggests a violation of causality, yet our solution is manifestly causal. How can this be? 
It turns out that these fluxes appear because of how the flow is divided into a mean and 
fluctuating part. Consider a simple example: an accretion flow consisting of a turbulent, 
unmagnetized fluid with negligible density fluctuations. At a given event, the flow has 
angular momentum I + 51 and radial four-velocity u'^ + Svl' . The mean flow is deflned so 
that {51) = and {5u^) — 0, where the brackets denote an average over t, 9, (j). Then the 
outward flux of angular momentum is 

T/ ^{p + u + p){l + 5l){u' + 5u'), (13) 

using the deflnition of the perfect fluid stress tensor. Averaging, 

T/ = {p + u+p) [lu' + {6l6u')] . (14) 

The flrst term in brackets is due to the mean flow; the second term is what we have called 
the "viscous" angular momentum flux. Evidently the outward flux of angular momentum 
(or energy; the same considerations apply to T^*") appears because of how we have divided 
the flow into a mean and fluctuating part. Correlations in the fluctuations merely bias the 
angular momentum of accreting fluid elements. Causality is preserved. 
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5.4. Assumptions and Limitations 

Finally, it is worth offering a frank; discussion of the assumptions behind our solutions 
and their limitations. A somewhat hidden assumption is that the accreting plasma 
is a two-temperature plasma with proton temperature much greater than the electron 
temperature. Only then will / ~ 1. This can be true only if most of the "viscous" 
dissipation goes into the protons and there is no collective effect that efficiently couples the 
protons and electrons. No such effect has yet been convincingly demonstrated to exist. 

The assumption that / is constant with radius is only likely to hold true if / ~ 1 
throughout the flow. We have nonetheless calculated solutions for constant / < 1 in order 
to illustrate the effects of a smaller value of / on the dynamical aspects of the flow. In 
the future, it is clear that models will need to include cooling processes and calculate / 
self-consistently. This task is complicated by the importance of Compton cooling in these 
flows, which depends not only on the local conditions at a particular radius, but also on 
the incident photon flux from all other radii. Thus far, cooling has been included in some 
detail in some models, but no model has included both a full treatment of cooling and fully 
relativistic dynamics. 

Some additional assumptions and limitations of our solutions are tied to our treatment 
of angular momentum transport. For example, we have vertically averaged the flow. This 
should produce a reasonably reliable solution close to the midplane, but is not predictive 
for flow near the poles. We could produce a full axisymmetric steady-state solution, but 
because it is not known how turbulent angular momentum transport varies with height 
in the accretion flow, such an effort would not significantly improve the reliability of the 
solution. In addition, the flow has been assumed smooth and steady. This is hkely to 
be true only in a time-averaged sense. In particular, if a ~ 1, then the turbTilcncc that 
transports angular momentum will be only marginally subsonic and shocks and substantial 
density variations are likely. 

6. Summary 

We have presented steady-state solutions in the Kerr metric for advection-dominated 
accretion onto black holes. The solutions extend from 2 x lO^GM/c^ down to just outside 
the event horizon at (1 -|- \/l — a'^)GM/c^. The flow passes through a sonic point and a 
viscous point which arises due to our use of a causal viscosity prescription. We find large 
variations in the character of the flow by varying the solution parameters, which include 
the black hole spin a, the viscosity parameter a, the advected fraction /, and the adiabatic 
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index 7. We find smooth solutions without shocks throughout the parameter space. 

The innermost part of the flow near the horizon becomes much hotter and denser as 
a approaches 1. Also, the amount of angular momentum transferred to the black hole, 
which is an eigenvalue of our solutions, drops as a increases. The structure of the flows 
is sensitive to a, the efficiency of angular momentum transport. For a < 10"^, the fiows 
exhibit a pressure maximum and an approximately constant specific angular momentum 
profile, similar to the "thick disk" or "ion tori" models. The same holds true for solutions 
with small values of /, even though the low-/ solutions arc geometrically thin and are 
similar to the classical thin disk solutions in a number of respects. 

Wc find that an equilibrium spin rate a^q exists for black holes accreting from an 
advcction-dominatcd flow. This results from the decrease in the angular momentum 
accretion rate as a increases. At aeq, the accreted angular momentum is just sufficient to 
counteract the effects of mass accretion and keep a constant. For large values of a and /, 
~ 0.8 — 0.9, well below the maximum value of a ~ 0.998 allowed by photon transport. 

We thank R. Narayan for helping to initiate this project and for his support, insight 
and encouragement, and H. Fak:kc, J. -P. Lasota, M. Rees, and M. Abramowicz for helpful 
discussions. This work was supported by grants NASA NAG 5-2837 and NSF AST 9423209. 
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Fig. 1. — Advection-dominated accretion flows for various values of a, aU with a — 0.1, 
/ = 1, and 7 = 1.4444. The solutions shown have a = -0.999, -0.9, -0.5, 0, 0.5, 0.9, and 
0.999. The six panels show the radial velocity V, density p, angular velocity Q, specific 
angular momentum Irj, temperature T, and vertical scale height as a fraction of radius H/R. 
On the panel showing V, the position of the sonic point is denoted by a filled circle, and 
that of the viscous point by a filled square. On the panel showing fl, the value of oo at the 
horizon is given by a dashed line. 

Fig. 2.— Similar to Fig. 1, but for flows with a = 0.001, 0.003, 0.01, 0.03, 0.1, and 0.3, aU 
with a = 0, / = 1, and 7 = 1.4444. 

Fig. 3.— Similar to Fig. 1, but for flows with / = 0.01, 0.03, 0.1, 0.3, and 1, aU with a = 0, 
a = 0.1, and 7 = 1.4444. 

Fig. 4.— Similar to Fig. 1, but for flows with 7 = 1.3333, 1.4444, 1.55, and 1.66, aU with 
a = 0, a = 0.1, and f — 1- 

Fig. 5. — (a) The angular momentum accretion rate j as a function of the black hole spin 
a, for a = 0.001, 0.003, 0.01, 0.03, 0.1, and 0.3, all with / = 1 and 7 = 1.4444. The dashed 
line denotes j = 2a; for solutions with j > 2a, accretion spins up the black hole to larger a. 
(b) same as (a), but for / = 0.01, 0.03, 0.1, 0.3, and 1, all with a = 0.1 and 7 = 1.4444. 

Fig. 6. — The equilibrium value Ugq, where j = 2a, as a function of a (top panels) and 
/ (bottom panels). The left panel in each case shows a^q, while the right panel shows 
log(l -aeg). 

Fig. 7. — The light lines show isodensity contours in the zu — rsinO, z — r cos 9 plane, where 
r,9 are the usual Boyer-Lindquist coordinates. Solutions are shown for a — 0.1, 10~^ and 
a = 0, 0.99. The heavy solid lines show the event horizon (inner curve) and the boundary of 
the ergosphere (outer curve). 

Fig. 8. — The Bernoulli parameter Be — rj£ — 1 as a function of radius for the four sets of 
solutions shown in Figs. 1-4. The four panels show solutions for various values of (a) a, (b) 
a, (c) /, and (d) 7. 

Fig. 9. — The outward viscous energy flux, — for solutions with a = —0.999, -0.9, -0.5, 
0, 0.5, 0.9, 0.999 and otherwise standard parameters. Notice that for a ^ 0.7 there is an 
outward viscous energy flux at the horizon. 
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